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Comment on " Electron-phonon Vertex in the 
Two-dimensional One-band Hubbard Model" 

In the past months a variety of experiments have 
pointed out an important role of the electron-phonon 
(el-ph) interaction in many physical properties of the 
cuprates 0, Q, 0] • These recent findings have triggered 
a renewed interest for a theoretical understanding of 
the electron-phonon properties in strongly correlated sys- 
tems. 

One of the most remarkable effects of the strong elec- 
tronic correlation on the el-ph properties is a to favor for- 
ward (small q) scattering in the electron-phonon vertex, 
q being the exchanged phonon momentum. This feature 
was investigated in the past by means of analytical tech- 
niques based on slave-bosons or Hubbard X-operators 

la, |(| . The assumption of forward scattering predomi- 
nance within an el-ph framework was shown to explain in 
a natural way several anomalous properties of cuprates 
as the difference between transport and superconducting 
el-ph coupling constants la, the linear temperature be- 
haviour of the resistivity J7j , the <i-wave symmetry or the 
superconducting gap [y, |8j • Small q scattering selection 
was shown moreover to be responsible in a natural way 
for high-T c superconductivity within the context of the 
nonadiabatic superconductivity |9j. 

So far, this important feature was analyzed only by 
means of analytical approaches in the U = oo limit and 
a definitive confirm of it based on numerical methods was 
lacking. With these motivations in a recent paper Huang 
et al have address this issue in the twodimcnsional Hub- 
bard model with generic U by using Quantum Monte 
Carlo (QMC) techniques on a 8 x 8 cluster [ljj. Their 
results provide a good agreement with the previous an- 
alytical studies and represent an important contribution 
to assess the relevance of el-ph interaction in correlated 
system. 

Moved by similar motivations, we have recently inves- 
tigated the onset of forward scattering predominance in 
the el-ph vertex by using slave-boson techniques based 
on four auxiliary fields |ll|. While the previous analyt- 
ical studies were limited to the U — oo limit 0, 0, @ , 
we were able in this way to evaluate the small momenta 
selection in the whole phase diagram of parameters U-n, 
where n is the electron filling, in direct comparison with 
Huang et al. |l0| . Our findings reproduce in more than 
a qualitative way the results of Ref. 01 1 an d they will 
be object of a future publication which is in preparation 
at the moment |l^ . 

Besides the direct numerical confirmation of small mo- 
menta predominance in correlated systems, Huang et al. 
report in their paper an increasing of the el-ph vertex as 
function of U in the large U regime for small phonon mo- 
menta and a relative saturation for large q |lQ] . The ori- 
gin of this puzzling feature has been not well understood 
so far and a possible connection with charge excitations 



driven by the exchange term J oc t 2 /U in the large U 
limit was mentioned [lOj. 

In this Comment we would like to make clearer the 
issue. In particular our analytical calculations suggest 
that the increase of the el-ph vertex function at small q as 
function of U is a peculiar feature of the high temperature 
regime and it could be associated to the tendency towards 
a phase-separation instability. 

Our technical approach was based on the four slave- 
boson method first introduced in Ref. . The electron- 
phonon vertex function for a finite U Hubbard model has 
been evaluated as linear response to an external field cou- 
pled with charge density. In the spirit of functional inte- 
gral representation the screening of the electron-phonon 
function is due to the Gaussian fluctuations of the auxil- 
iary boson fields around the mean-field solution. A simi- 
lar study at finite U was introduced by A. Lavagna who 
however did not addressed the momentum modulation 
of the electron-charge density response ^ij. A different 
choice for the renormalization of the ^-operators was in 
addition done in comparison with Refs. |l3j, [lj] to avoid 
unphysical divergences in the electron-slave boson ma- 
trix elements. Technical details will appear in Ref. [l2| . 
Based on this analytical approach, we can now evaluate 
the electron-phonon function with the same physical pa- 
rameters of Ref. ^jj , namely n = 0.88, (3 = 2, where n is 
the electron filling and (3 = t/ksT the inverse tempera- 
ture in unit of the nearest neighbor hopping parameter t. 
As only slight difference with respect to Ref. [ljj we as- 
sume the incoming electron momentum p to be averaged 
over the Fermi surface and the exchanged frequency to 
be exactly zero. These marginal differences are expected 
to not significantly affect the comparison between our 
results and Ref. |l0| . 

In Fig. we plot the electron-phonon vertex function 
<7(q) at q = (7r/4, tt/4) and q = (tt, 7t) as function of the 
Hubbard repulsion U for n — 0.88 and (3 = 2. For q = 
(tt/4, 7r/4) we note that while for relatively small U the 
el-ph vertex function is steady decreasing with U, such 
a behaviour has an upturn for U ~ 8 until a divergence 
occurs for U c ~ 9.3. According this view one is attempted 
to associate the upturn of <?(q) as function of U as an 
incipient transition towards some charge instability. Note 
also that for q = (tt, tt) no charge instability is observed. 

The appearance of charge instability can be also de- 
tected by looking at g(q) plotted as function of q (q = 
(q,q)iT). The evolution of g(q) by varying the Hubbard 
repulsion U is shown in Fig. H^c where we see that the 
el-ph vertex function is initially suppressed at small q by 
increasing U (panel b) , then it increases as function of U 
(panel c) until a divergence is established. Note that for 
U = 9 a lattice instability already occurs although for a 
q (~ 0.2) less than q ~ 1/4. This reflects the fact that 
a instability for q = (phase separation) is first onset 
for some critical value of U, and then the vector insta- 
bility is gradually shifted by further increasing of U. In 
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FIG. 1: (a): Electron-phonon vertex <?(q) as function of the 
Hubbard repulsion U for q = (7r/4, 7r/4) and q = (7r,7r). The 
temperature was set here f3 = 2 and the electronic filling 
n = 0.88. (b)-(c) Plot of <?(q) as function of q = {q,q)ir 
for U = 2,4,6 (respectively circles, squares and diamonds in 
panel b) and for U = 7, 8, 9 (respectively crosses, triangles up 
and triangles down in panel c). 
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FIG. 2: (a): Electron-phonon vertex <?(q) at q = (7r/4,7r/4) 
as function of the Hubbard repulsion U for (5 = 2 and (3 = 50. 
Electronic filling n = 0.88. (b) Phase diagram for phase sep- 
aration (q = 0) and for the charge ordering q = (tt/4, 7t/4) 
instability in the T-U space. The dashed line represent the 
band narrowing factor Z corresponding to the phase separa- 
tion instability line. 



this perspective it is not surprising that momenta on the 
Brillouin zone edge q = (tt, 7t) are less sensitive to the 
increase of U. 

The similarity between our findings and the Quantum 
Monte Carlo analysis suggests that also the upturn of 
<7(q) as function of U reported in Ref. could be re- 
lated to the same tendency towards phase separation or 
charge instabilities. This does not imply however that 
phase separation is effectively established, and it should 
be remarked that the actual occurrence of phase separa- 
tion in the Hubbard model is still object of debate 
On one hand expansions around the mean-field solution 
even including Gaussian fluctuations could enforce un- 
physical instabilities which could disappear once higher 
order fluctuations are taken into account, especially in 
two dimensional systems. On the other hand, small size 
cluster effects (L x L) and large temperature effects in 
Ref. ^3 question in principle the generalization of the 
QMC results in the thermodynamic limit (L — > oo) and 
at low temperatures. Our results should thus viewed as 
analytical indications which can trigger further numerical 
work. 

As a final step, we can also employ our slave-boson 



analysis to extend the range of inve stig ation in regions of 
parameters not addressed in Ref. [lfj. In particular we 
show that a crucial role is played by the temperature that 
is limited in QMC techniques by the sign problem and 
by the requirement to be larger that the energy spectrum 
discretization. 

In Fig. we report the dependence of the el-ph vertex 
function g(q) at q = (it/4, tt/4) as varying U in the large 
(f3 = 2) and small ((3 — 50) temperature limit. As a 
surprising result, no charge instability is found in the 
low temperature regime in the contrast with the high 
temperature range ((3 = 2). Same conclusions hold true 
for the phase separation q = instability (not shown in 
figure) . This results in thus compatible with the absence 
of phase separation at zero or low temperature. 

In order to understand in more detail the origin of the 
phase separation instability as function of the tempera- 
ture T we show in Fig. [2b the phase diagram in the T vs. 
U space with respect to phase separation (q = 0) and to 
charge ordering (q = (7r/4, 7t/4)). As above mentioned, 
a phase separation instability occurs in our slave-boson 
calculations only above a certain temperature T/t ^ 0.2 
((3 % 5) 0|- As expected, finite q instability, in the ab- 
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sence of any long-range Coulomb repulsion, is prevented 
by the occurrence of phase separation at q = in the 
whole phase space. A interesting insight comes from 
the comparison of the critical temperature T c at which 
the instability towards phase separation occurs with the 
band narrowing factor Z due to the correlation effects 
[H E (dashed line in Fig. 0d). The similar de- 
pendence on U of T c and Z points out that the onset of 
phase separation by increasing temperature is ruled by 
the comparison between the temperature T and the ef- 
fective bandwidth W = Z8t energy scales. In particular 
phase separation is established when T become is of the 
same order of with Zt (= W/8). Once again we stress 
that the phase separation instability found by our slave- 
boson calculations which include Gaussian fluctuations 
could be washed out when higher order fluctuations are 
taken into account, so that it should be regarded as in- 
dicative of tendency towards this kind of instability. Nu- 
merical work based on Quantum Monte Carlo techniques 
will help to answer about the effective role of phase sep- 
aration or charge ordering in the Hubbard model. 
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